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Abstract 

We consider infinite random casual Lorentzian triangulations emerg- 
ing in quantum gravity for critical values of parameters. With each 
vertex of the triangulation we associate a Hilbert space representing a 
bosonic particle moving in accordance with standard laws of Quantum 
Mechanics. The particles interact via two-body potentials decaying 
with the graph distance. A Mermin-Wagner type theorem is proven 
for infinite-volume reduced density matrices related to solutions to 
DLR equations in the Feynman-Kac (FK) representation. 
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1 Introduction 



In this paper we prove a Mermin- Wagner (MW) type theorem (cf. [15] ) 
for a system of quantum bosonic particles on an (infinite) random graph 
represented by a causal dynamical Lorentzian triangulation (in brief: CDLT). 
The CDLTs arise naturally when physicists attempt to define a fundamental 
path integral in quantum gravity. The reader is referred to p2] for a review of 
the related publications and to [H] for a rigorous mathematical background 
behind the model of CDLTs. More precisely, we analyze a quantum system 
on a random 2D graph T generated by a natural "uniform" measure on the 
CDLTs corresponding to a "critical" regime (see below). 

In a modern language, the spirit of the quantum MW theorem is that in a 
2D lattice model (more generally, for a model on a countable bi-dimensional 
graph), any infinite- volume Gibbs state (regardless of whether it is unique 
or not) is invariant under the action of a Lie group G provided that the in- 
gredients of local Hamiltonians are G-invariant; see [15J. These ingredients 
include the kinetic energy part, the external field and the interaction po- 
tential. The mathematical constructions used for the proof of this theorem 
require a certain control over these ingredients: compactness of a configura- 
tional space associated with a single vertex of the lattice (or a more general 
bi-dimensional graph), smoothness of the interaction potential (or its essen- 
tial part), decay of the interaction potential for large distances on the lattice 
(or on the graph), "regularity" of the lattice (graph) geometry. 

A principal question that needs a careful consideration is about the def- 
inition of a quantum Gibbs state in an infinite volume. For the so-called 
quantum spin systems, with a finite-dimensional phase space of a single spin 
(and consequently, with bounded local Hamiltonians), such a definition is 
given within the theory of KMS (Kubo-Martin-Schwinger) states; see, e.g., 
[1]. A version of the MW theorem for a model of this type on a 2D square 
lattice had been established in [I], [16] and subsequent publications. The 
KMS-based results (under suitable aforementioned assumptions) can be ex- 
tended to the model of classical spins on a random graph of the type consid- 
ered in the present paper; cf. pU]. However, the KMS-theory is not efficient 
for the case of interacting quantum particles where the one-particle kinetic 
energy operator is equal to —A/2 (A stands for a Laplacian on a compact 
manifold). This is a standard quantum-mechanical model, and the fact that 
the concept of infinite-volume Gibbs state did not receive so far a properly 
working definition for such a system was perceived as a regrettable hindrance. 
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In this paper we adopt the definition of an infinite-volume Gibbs state 
(more precisely, of an infinite-volume reduced density matrix (RDM, for 
short)) from the papers [8], [9]. In these papers, a class of the so-called 
FK-DLR states has been introduced, and an MW theorem was established 
for quantum systems on a bi-dimensional graph T ~ (V, £) where V = V(T) 
is the set of vertices and £ = £{T) the set of edges. As was said above, 
in the present paper we deal with a random graph T (a CDLT for critical 
values of parameters). After checking that a typical realization T = of 
the random CDLT satifies required properties, we use the constructions from 
[H], [H] (going back to [I], [TU]) and prove the main results of the paper (see 



Theorems 2.1, 2.2 and Theorems 3.1, 3.2). 
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2 Basic definitions 



2.1. Lorentzian triangulations in a critical phase. 

The graphs under consideration is a triangulation T of a cylinder C = 
S x [0,oo) with the base S, a unit circle in M 2 . Physically, C represents a 
(1 + l)-dimensional space-time complex. Geometrically, C can be visualized 
as a complex plane C with a family of concentric circles {z : \z\ = n}\ 
here the origin z = is treated as a "circle" of infinitesimal radius. The 
following properties of T are assumed: each triangle belongs to some strip 
S x [£,£ + 1], £ — 0, 1, 2, . . . such that either (i) two vertices lie on S x {£} 
and one on S x {£ + 1} or (ii) two vertices lie on S x {£ + 1} and one on 
5 x {£} and exactly one edge of triangle is an arc of a circle S x {£} in the 
case (i), or S x {£ + 1} in the case (ii). In case (i) we speak of an upward 
triangle, or simply up-triangle, and in case (ii) of an downward triangle, or 
down-triangle. This includes also a "degenerate" picture where two vertices 
of a triangle coincide, and corresponding edge forms the circle. For £ = it is 
requested that the graph under consideration generates a degenerate picture 
(i.e., the graph has a single up-triangle in the strip S x [0, 1], see Figure [Ija). 
This particular triangle is called the root triangle, and its side represented 
by the edge cutting through the strip is called the root edge. Moreover, the 
(double) vertex lying on S x {0} is called the root vertex. 

Finally, we consider graphs modulo an equivalence (that is, up to a home- 
omorphism of C preserving all circles S x {£}, £ = 0, 1, 2, . . .). 

Definition 1 A rooted infinite CDLT is defined as an equivalence class of 
(countable) graphs with the listed properties, under the above equivalence. 
Depending on the context, we use the notation T = for a representative 
or the whole equivalence class of graphs involved, and speak of the vertex set 
V = V(T) and the edge set £ = S(T) in the same fashion. 

A similar definition can be used to introduce a rooted CDLT on a cylinder 
= S x [0,N] (a rooted CDLT of height N). The corresponding notation 
is T N ~ (V(T N ),S(T N )) or even T N ~ (Vn,£n)- As in [TOj, we denote by 
LTjv and LT^ the sets of CDLTs on Cjy and C respectively. 

To introduce a probability distribution on hT^ and ultimately on LToo, 
we use a special 1-1 correspondence between the rooted CDLTs and rooted 
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Figure 1: (a) An example of Lorentzian triangulation. Some (not all) up- 
and down-triangles are marked, (b) The triangulation T is parametrized by 
the spinning tree T, which is represented by bold lines. 

trees (that is, graphs without cycles and with distinguished vertices). Namely, 
we extract a subgraph in T by selecting, for each vertex v G T, the leftmost 
edge going from v downwards and discarding all other edges going from v 
horizontally or downwards, see Figure |T|(b) . The graph T C T thus obtained 
is a spanning tree of T, cf. [3] and [TJ]. Moreover, if one indicates, for each 
vertex of T, its height in T then T can be completely reconstructed when we 
know T. We call the correspondence T -H- T the tree parametrization of the 
CDLT. It determines a one-to-one bijection m between the set LT^ and the 
set of infinite rooted trees Too'- 

We will use the same symbol m for the bijection TJv -H- LT^r where 7jv is the 
set of all planted rooted trees of height N . 

By virtue of the tree-parametrization, we will specify a probability dis- 
tribution on CDLTs by specifying a distribution defined on trees. More 
precisely, suppose P^ ee is a probability measure on Tn- Then the measure 
on LTtv is determined by 

pL T (m(T)) = P£ ee (T), V T G Tn- 

Conversely, let P tree be a probability distribution on Too- Then the distribu- 
tion P LT on LTqo is given by 

P LT (m(A)) = P tree (A), V A G T(%o), 
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where J 7 (Too) is standard cr-algebra generated by cylinder sets. In future we 
omit indices in the notation for the distribution P. 

To constract a critical CDLT model we define the corresponding measure 
on Too related to a critical Galton- Watson (GW) process £. For this aim, 
we set fi = {pk} to be an offspring distribution on N = {0, 1, . . . } with 
mean 1, and define the critical Galton- Watson (GW) branching process £. 
Conditional on the event of non-extinction, the GW process becomes the 
so-called size-biased (SB) process £ = {k n , n = 0, 1, 2, . . .}. In our context, 
k n yields the number of vertices on the circle S x {n} in the random infinite 
rooted CDLT. The reader can consult [T5] for the formal background for SB 
branching processes. 

In particular, the distribution of an SB process £ is concentrated on the 
subset S of To formed by the so-called single-spine trees. A single-spine tree 
consists of a single infinite linear chain so,s\, . . . called the spine, to each 
vertex Sj of which there is attached a finite random tree with their root at 
Sj. (Here Sq is the root vertex of the whole tree.) Furthermore, the gener- 
ating function for the branching number v at each vertex Sj is f'(x) where 
f(x) is the generating function of the initial offspring distribution //. More- 
over, the individual branches are independently and identically distributed 
in accordance with the original critical GW process |13j . 

Let a 2 stand for the variance of the offspring distribution v. Then 

E(k n \k n _ 1 ) = k n . 1 + a 2 . (2.1.1) 

In fact, let v = {pk} be the SB offspring distribution with p k = kp k . Then the 
distribution for k n conditioned upon the value k n _i is identified as follows (cf. 
|13j). We choose at random one particle among particles and generate 
the number of its descendants according to the distribution v with mean 
a 2 + 1. According to (2.1.1), the number of descendants for other k n _\ — 1 
particles is generated, independently, by the distribution \i. 

Throughout the paper we assume that the offspring distribution \i has 
the mean 1 with finite second moment. Let P be the corresponding SB 
Galton- Watson tree distribution. 

2.2. The local quantum Hamiltonians on CDLTs. Let M = R d /Z d 

be a unit (^-dimensional torus with the flat metric and the induced volume v. 
A basic quantum model uses the Hilbert space "H = L 2 (M, v) as the phase 
space of a single quantum particle. The single-particle Hamiltonian H acts 
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in TL as the sum: 



(H(f))(x) = --(A(f))(x) + U(x)(f)(x), x G M, (j>en. 



(2.2.1) 



Here A is the Laplacian on M and the function U : x e M i— >■ R gives an 
external potential. Under the assumptions (2.3.1) upon U adopted in this 
paper, if is a self- adjoint operator bounded from below and with a discrete 
spectrum such that V (3 > 0, exp [—f3H] is a (positive definite) trace class 
operator. 

Given a CDLT T we use the notation Tjy for the subgraph in T with 
the set of vertices V N = V(T N ) of the form V N = V(T) n{Sx{0,... N}} 
and the set of edges £(Tn) = £(T) fl (Vn x Vat). The phase space of a 
(bosonic) quantum system in Vn is the tensor product "Hat = "H® ViV ; an 
element 0a? G "Hat is a function 



x(7V) = x ViV = {a;(i),i e V N } e M xVn ^ C 



(2.2.2) 



square-integrable in dx(iV) = [ [ i;(dx(i)). The Cartesian power M xVjv can 

be considered as the configuration space for the classical prototype of the 
quantum system in V^. 

The local Hamiltonian of the system in Vn acts on functions <fi e 
n® VN : given x(iV) = (x(j), j e V N ) G M xVn , 



(H N <j>)(x(N)) 



W,f))V(x(j),x(f)) 

.*eV]v j,j'eV N xV N 



0(x(iV)), 
(2.2.3) 

where J(d(j,j')) V(x(j),x(j')) represents the interaction between spins x(j) 
and x(f) at sites j and j'. Next, H(i) stands for the copy of operator H 
acting on variable x(i) e M and d(j,f) for the graph distance from vertex j 
to f. 



A more general concept is a Hamiltonian H 



N\x c (N) 



in the external field 



generated by an (infinite) configuration x c (iV) = {x(j'), j' e Vn} £ M xVjv 
where Vn = V(T) \ Vat. As before, operator H n \x c (n) acts in "Hat: given 
(f)en N and x(JV) = (x(j), j e V N ) G M xV ^, 



(^1x^(^)0) (x(iV)) 



^+ ^ J(d(j,j'))V(x(j),x(f)) 
(j,j')ev N xV N 



0(x(A)). 
(2.2.4) 
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Again, under assumptions (2.3.1)-(2.3.2) upon J and V, H N and Hn\ x c (n) 
are self-adjoint operator bounded from below and with a discrete spectrum 
such that V (3 > 0, Gp t N = exp [—/3Hn\ and Gp t N\^(N) = exp [ — (3H n \^(n)] 
are (positive definite) trace class operators. 

The operators Gp^ and Gp^\^(N) are called the Gibbs operators (in 
volume Vat for the inverse temperature (3 and, in the case of Gp ; N\^^N), with 
the boundary condition x c (iV) = Xy ). The traces 

^/3,n = tTy. N Gp t N and z-p,N\x?(N) = toy. N Gp t N\x?(N) (2.2.5) 
give the corresponding partition functions. The normalized operators 

Rp,N — ^ — Gp tN and Rp t N\x. c (N) = ^ Gp tN \^( N ) (2.2.6) 

^P,N ^p,N\x c (N) 

are called the density matrices (for the corresponding Gibbs ensembles); these 
are positive definite operators of trace I. Given n G {0, . . . , N}, the partial 
traces 

r< £n = tr n N \ n Rp,N and Rp n ^ NlxC{N) = ton N \ n R/3,N\x°(N) (2.2.7) 

yield positive definite operators Rp n \ and Rp n ^ N \^(N) m "^«> °f trace 1. Here 
Hn\u stands for the Hilbert space %®( v jv\ v «). These operators are called 
the reduced density matrices (RDMs). Note the compatibility relation: V 
< n < ri < N: 

R p n ]v = tT n n ,\ n Rp n t N and i?^| xC(JV ) = tr 7V\»^&v|x«(A0 (2.2.8) 



2.3. Assumptions on the potentials. The group of symmetries. 

We suppose that the potential U has continuous derivatives whereas V has 
continuous first and second derivatives: V Ob ■ Ob ■ 0C G M 

\U(x)\,\V x U(x)\<U, (2.3.1) 

\V(x',x")\, \V x ,V{x',x")\ , \V x „V{x',x")\ , \V x ,V x „V{x',x")\ < V, (2.3.2) 

where U, V G (0, oo) are constants. 

Next, suppose that ad'xd matrix A is given, of the row rank d! where d' < 
d. We consider a cf -dimensional linear group G acting on M and preserving 
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the volume v: (g, x) G G x M i— > gx G M. More precisely, g is identified with 
a real d'-dimensional vector 9 — ... ,9^) and the action is given by 

gx = x + 6A modi. (2.3.3) 

Remark 2.1 The group G can fre compact (in which case G is a torus of 
dimension d! ) or non-compact (then G is M. d ). 

We assume that the functions U (x) and V(x, x') are invariant with respect 
to the group G: Vg G G and x,x' G M 

U(gx) = U(x), V(gx,gx') = V(x,x'). (2.3.4) 

Finally, we assume that the function r G (0, oo) h- > J(r) in (2.2.3) and 
(2.2.4) is a bounded monotone decreasing function satisfying the condition 

j{r) -(vL)^ r - 2 - (2 - 3 - 5) 

These assumptions are in place throughout the paper. 

As usually, the action of the group G generates unitary operators in H: 

S(g)(f)(x) = (pig^x), xeM, (f>en. (2.3.6) 

Let S( N \g) be the tensor power of S(g) which acts in 1-Ln'- for any <pN G Hn 

S {N) (g)M*(N)) = Mg-'AN)), (2.3.7) 

where x(iV) = {x(i), i G V N } G M xVn and g _1 x(iV) = {g^x(i), i G V^}. 

2.4. Limiting RDMs in an infinite volume. We are interested in 
the 'thermodynamic' limit iV — > oo. Ideally, one would like to establish a 
convergence of the RDMs Rp% and ^ n L x c (Jv) ^° a nm hing RDM in "H n 
as — > oo. A suitable form of convergence is in the trace norm in T-t n , 
guaranteeing that the limiting operator is positive-definite and has trace 1. 
A more general question is whether the families {R^ N } and {R^ N ^c^} are 

compact. If we manage to check that {R^ N } and {R^ N \ x cr N ^} are compact 
families for any given n then, invoking a diagonal process, we can consider 
a family of limiting RDMs {R^, n = 0, 1,2, . . .} (in the case of operators 
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Rp h\x c (N) ^ ne lifting RDMs may depend on the choice of the boundary 
conditions x c (iV)). The consistency property (2.2.8) will be inherited in the 
limit: V < n < n' < N, 

A consistent family of RDMs i?L defines a state of (i.e., a linear positive 
normalized functional on) the quasilocal C*-algebra constructed as the clo- 
sure of the inductive limit of 53 n as iV — > oo where 23 at is the C*-algebra of 
the bounded operators in T-L N , cf. pQ. This motivates a study of properties of 
limiting RDM families {R^}. Our results in this direction are summarised 
in Theorems 12.11 and 12.21 



Theorem 2.1 Fix (3 > 0. For P-a.a. CDLT T 6 %o, Vn = 0, 1,2,..., the 
family of the RDMs {Rp N , N = 1, 2, . . .} is compact in the trace norm in 

7i n . Similarly, {Rp N i x cnv) ; N = 1, 2, . . .} is a compact family V choice of 
the boundary conditions x c (iV). 



Theorem 2.2 Let Rf } be any limiting-point operator for the family {Rp N \ xC( 
iV = 1, 2, . . .}. Then, V g 6 G, operator S^ n \g) commutes with Rp : 

Rf ) = S^(g)R i ; ) (S^(g)y 1 . (2.4.2) 



Remark 2.2 The statement of Theorem 2.2\ is straightforward for the limit 



points Rp of the family {Rp N , N = 1, 2, . . .} but requires a proof for the 
family {R^ N]jc c {N) , 



iV = l,2,...}. 



The main role in the proof of Theorems 2.1 and 2.2 is played by the 



Feynman-Kac (FK) representation for the RDMs Rpjj and Rp riJ m x c/ N \ anc ^ 

their limiting counterparts R^ . This representation is discussed in the next 
section. 



(nj 



10 



3 The FK ensembles of paths and loops 



3.1. The FK representation for the Gibbs operators. The Gibbs 
operators G^,jv an d Gp t N\ x c (N) ac t as integral operators, with kernels Kp^ 

and Kp t N\x c (N) 

(G? AJV ^)(x(JV))= J Kf,, N (x{N),y{N))My(N))dy{N), 

M xV N 

(G AJV|x c (JV) ^)(x(JV)) = y ^| x c (JV) (x(iV),y(iV))0 iV (y(iV))dy(iV). 

M xV iV 

(3.1.1) 

Here y(iV) = {?/(«) : i G Vat} and we use a shorthand notation dy(iV) = 

n v(Mi))- 

ieVjv 

Further, the kernels Kp tN and Kp iN \^(N) admit the FK-representations 



summarized in Lemma 3.1 The proof of this lemma follows the standard 
lines and is omitted. The reader can confer [6] for details. 

Given points x, y G M, let W^ y denote the space of continuous paths 
u = UJ X) y : r G [0, (3] i— > uJ(r) G M, of time-length /?, beginning at x and ter- 
minating at y. Next, let Pf stand for the (unnormalized) Wiener measure 
on VKf„, with ^(W^) = p l3 (x,y) where pP(x,y) is the value of the transi- 
tion density from x to y in time /3. Furthermore, given particle configurations 
x(iV) = {x(i)},y(N) = {y(i)} G M xVn , we set: 

In other words, an element U(N) = oJ*.(N),y{N) G W^x(2V)y(JV) * s re P resen ted 
by a collection of paths {co x (i),y(i}} where ^x(i),y(i) e We ca ^ sucn a 

collection a path configuration over Vat. Moreover, under measure Pwjv) yW) 
the paths u^i),^) are independent and each of them follows its own marginal 
measure P ,.s ,. v 

Further, we need to introduce functionals h(cJ(N)) and fo(u>(iV)| x. c (N)) 
describing an integral energy of the path configuration u(N) and its energy 
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in the potential field generated by x c (iV): 



h(p(N))= h^(uJii),uj(i')) (3.1.3) 

(i,i')<=V N xV N 



Namely, for i — i' and u> G 



where h %,% ' \u(i) : uj(i')) represents an integral along trajectories u(i) and uj(i'). 

x(i),y(i)' 

h^(cd,cd) = J drU(uj(r)). (3.1.4) 



o 

and for i ^ %' and to G Wl? M „ M , uj' G Wt- n , 



p 

h^{u,u') = J(d(i,i')) J drV(cu(r),cJ'(T)). (3.1.5) 

o 

Pictorially, h l ' l {u(i),uj(i)) yields an energy of path cJ(i) in the external field 
gernerated by the potential U and /i M (uj(i),UJ(i f )) the energy of interaction 
between paths uj{i) and w(i'). Accordingly, h(u(N)) gives a full potential 
energy of the path configuration cJ(N). 
Similarly, 

h(UJ(N)\x c (N)) = h(p(N)) + h^\u{i)^ c {i')) (3.1.6) 

iev N ,i'ev N 

13 

where h^'(uJ{i), x c (i')) = J(d(i,i'))JdrV(uJ(i,r),x. c (i')). 

o 

Lemma 3.1 The integral kernels if^jv^iV), y(iV)) and 
Kp jN i x c( N) (x(N),y(N)) are given by: 

K PtN ( X (N),y(N)) = J Pf W;y(iV) (du;(7V))exp [- h(u(N))] (3.1.7) 



VV x(N),y(N) 



and 



w p 



(3.1.8) 
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3.2. The FK representation for the partition functions and 
RDMs. Lemma [3. 1| implies a working representation for the partition func- 
tions H^jv and 'Bp ) N\xP(N) ( see (2.2.5)). More precisely, a key ingredient in 
the corresponding formulas will be the space W£ x = W@ of closed paths 
(starting and ending up at the same marked point x G M); we will employ 
the term a loop to make a distinction with a general case. Accordingly, the 
notation Pf x = Pf will be in place here. Note that measure P^ in essence 
does not depend on the choice of the point x G M. Furthermore, the nota- 
tion co = uj x G W% will be used for a loop with the marked initial/end point 
x, omitting the bar in the previous symbol u. Next, we set: 

W p im = x Wt„ P^n = x . (3.2.1) 

An element u>(N) G WffAO * s represented by a collection of loops {uj x (i)} 
where oj x ^ G W^y, such a collection is called a loop configuration over Vat. 
As before, under measure P^tv) * ne 1°°P S u x(i) are independent and each of 
them follows its own marginal measure Pf/. v 

At this point we apply the Mercer theorem guaranteeing that the traces 
tr nN G(3 tN and S /3iA r| x c (A r ) = tr Wjv G /3iA r| x c (Ar) are given by the integrals of the 
corresponding kernels Kp t isr(x(N),y(N)) and K^ ! jsr\^c^(-x.(N),y(N)) along 
the diagonal x(iV) = y(N). This leads to Lemma 3.2 below. 

Let us denote: 

/ du{N) := / dx(iV) f Pf (N) (du(N)). (3.2.2) 
J Jm xV n J 



Lemma 3.2 The partition functions H^jy an< ^ — /8,iV| x°(at) « r e gwen 6y: 

= y duj(N) exp [ - /i(w(JV))] (3.2.3) 



and 

^/3,N\x c (N) 



J du;(N) exp [ - h(u>(N)\ x c (iV))] . (3.2.4) 
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Let us now turn to the RDMs R^ N and Rp mxcrm- These operators are 
again given by their integral kernels: 

(fl£U0(x(n)) = j F^(x(n),y(n))^(y(n))dy(n), 

M xV n 

{ R th\x.*(N)<f>n)(x(n)) = J F^ |xC(7V) (x(n),y(n))0„(y(n))dy(n). 

M x v " 

(3.2.5) 

Next, Fg^-(x(n), y(n)) and F^-| xC ^(x(n),y(n)) are called reduced density 
matrix kernels (RDMKs). They can be written in the form 

F /3,7v( x H ; yH) = „ > 

— /9,iV 

^/3,JV|x c (iV) 

where quantities S^ r (x(n),y(n)) and H^| xC(Ar) (x(n), y(n)) admit 



(3.2.6) 



representations similar to (3.2.3) and (3.2.4), see Lemma 3.3 
We will use a notation similar to (3.2.2): 



du{N\n):= dx(iV\n) / P x(JVW) (du;(iV \ n)) (3.2.7) 



W P 

,V x(N\n) 



where x(iV \ n) stands for a particle configuration {x(j),j G \ V n } and 
u>(N\n) for the loop configuration {u;(j), j G Vjv\V n }. Symbol V will be used 
for concatenation of particle configurations and for concatenation of path and 
loop configurations (originally defined over disjoint sets). Accordingly, for a 
path configuration uJ(n) G Wf(n) y(n) over an< ^ a ^°°P configuration u(iV\ 
n) over VAr\V n , the energies h(<jJ(n)Vu}(N\n)) and /i(cJ(n)Vu>(iV\n)| x c (iV)) 
are defined as in (3.1.2)-(3.1.5). 

Lemma 3.3 The numerators sl"^-(x(n), y(n)) and 
S A^|xc(7V)( x ( n )>y( n )) are #™ en &?/; 

sW(x(n),y(n))= J Pf (n) , y(n) (da;(n)) 



x(n),y(n) 



(3.2.8) 



x y du>(iV \ n) exp 
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and 

^|WxW,yH) = / Pf {n)>y( „)(d^(n)) 

<n),y(n) (3-2-9) 

x y du(iV\n)exp [ - /i(w(n) V w(JV \ n)\ x c (iV))] . 



The proof of Lemma |3.3| consists in translating the partial traces into 
the integrals of the kernels and Fpm^crm of the operators R^ N and 

Rpm x c(xy We omit it from the paper. 

3.3. The FK-DLR equations. The representations (3.2.3)-(3.2.4) 
suggest introducing probability distributions /xjv and Hn\*. c {n) on loop con- 
figurations a; A?; with the densities (the Radon-Nikodym derivatives) 

Md^AQ) exp [-h(u,(N)j\ 

Pn{u{1\)) .— — - — - 

du(N) /ooiN 
, rAm ^| x c W (du;(Ar)) exp [-%,(iV)|x c (iV))] 
Pn\x?(n)(u{N)) := - 



dw(iV) 



-/3,AT|x c (Af) 



A crucial property is that the measures /i/v and /xjv| x c (jv) satisfy DLR (Dobrushin- 
Lanford-Ruelle)-type equations. Namely, let p^\oj{n)\(jj(N \ n)) and 
p$ xC , N J(jj(n)\u)(N \n)) stand for the conditional densities generated by 
and yUAr| x c(/v), respectively, for the loop configuration u?(n) over V n given a 
loop configuration uj(N \ n) over V/v \ V n . Then 

n Wu^LfAT\ „U - M^MM^X")) _ exp [ - h(u(n)\u(N \ n))] 

(n) / / \ I » fi N \xc(N)(fa>(n)\uj(N \n)) 
p { N ; xHn) (u(n)\u;(N\n)) := 

exp [ - /i(w(n)| w(iV \ n) V x c (iV))] 

S/3,n|x-(iV)(^(iV \ n)) 

(3.3.2) 

Here h(u)(n)\u)(N\n)) and /i(u;(n)| w(iV\n)Vx c (iV)) stand for 'conditional' 
energies and 'Ep tn (u)(N\n)) and 5 /f 3 )T1 | x c(jv)(u?(iV\n)) for 'conditional' parti- 
tion functions: 

h(u(n)\u(N \ n)) = h(u(n) V lj(N \ n)) - h(w(N \ n)), (3.3.3) 
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h(u(n)\u(N \ n) V x c (iV)) = h(u(n) V u(N \ n)\x. c (N)) 

-h(uj(N\n)\x c (N)), [ } 



E,3 >n (u(N \n)) = J du;(n) exp [ - h(ut(n))\w(N \ n)] (3.3.5) 

and 

Z MxHN) {u{N\n)) = J d w (n)exp [-h(u>(n)\ u{N\n) Vx c (iV))] . (3.3.6) 

We call Eqn (3.3.2) the FK-DLR equation in volume Vat. 
Concluding this section, we give an expression for the kernels Fg^(x(n), y (n)) 
and F^ |xC(JV) (x(n),yH): V < n < n' < N: 

FW(x(n),y(n)) = J d ^ N \ n ') €^^^ / M*'\n) 
x ! ^ Pf (n))y(n) (dw(n)) exp [-h(u(n) V u(n' \ n)\u{N \ n'))], 

x(n),y(n) 

Fp,N\ x c (iv)( x ( n )' y( n )) = 7 cL,(iV \ n \ w/) v xC(jV 



(*0) 

x exp [-/i(w(n) V u>(n' \ n)\ u(N \ n') V x c (iV))]. 



x / d«(n'\n) / ^ (n) , y ( B) (d5>(n)) 



(3.3.7) 

For n = n', the integral ^ du^n' \ n) is omitted. 

Our next goal is to write down FK-DLR equations for the whole of 

V = V(T). Here we consider a probability measure \i — \iy on infinite 
loop configurations f2 = fiy over V (for a formal background, see jH])- The 
equation is written for p( n \u)(n)\Q c (n)), the conditional probability density 
for a loop configuration u)(n) over V n , given a loop configuration f2 c (n) over 

V \ V„. This density should be given by 

pCVMIsrM) := = «p[ : M«Win-W)] (3 . 3 . 8) 

du>(n) ^,n(" (n)) 

Like /i (u> (n) | u? (N\n)) and Sg >n (u;(./V\7z)) before, the quantities h(uj(n)\Q c (n)) 
and S ( g i „(r2 c (n)) represent the conditional energy and the conditional parti- 
tion function. They can be defined as the limits 

h(u>(n)\fL c (n)) = Urn h(u(n)\n(N\n)), (3.3.9) 
16 



S An (n c (n)) = lim E p , n (Q(N\n)) (3.3.10) 

N— >oo 

where Q(N \ n) stands for the restriction of f2 c (n) to Vat \ V n . The existence 
of the limit will be guaranteed by the assumption (2.3.5) V u?(n) and f2 c (n) 
for P-a.a. T G Too- 

Formulas (3.3.7) admit a generalization to the infinite-volume situation: 
Vfl<n<n', 



x / Pf (n) , y(n) (d^H) exp [-h(U(n) V W (n' \ n)\ ft >'))]; 



w , 



(3.3.11) 

owing to the FK-DLR propety, the RHS in (3.3.11) does not depend on the 
choice of n' > n. Moreover, the integral 



J dx(n)Fj n) (x(n),x(n))=KV) 



Consider the operator in T-L{n) = %® Vn with the integral kernel 
Fg (x(n),y(n)) given by (3.3.4). The aforementioned properties imply that 
the trace tr^^i?!* 1 = 1 and the following compatibility relation holds true: 

Rp = tr W ( n /\ n )i?^ n . (3.3.12) 

Thus, were the operators positive definite, we could speak of an infinite- 
volume state of the quasilocal C*-algebra 23. Cf. Remark 2.2. Notwithstand- 
ing, we state our main result: 

Theorem 3.1 Under the above assumptions, any limit-point operator Ry^ 
from Theorem 2.1 is a positive definite trace-class integral operator of trace 1 
and with the kernel F^ admitting the representation (3.3.4 ) where probability 
distribution fi satisfies the infinite-volume FK-DLR equations (3.3.8). 

Theorem 3.2 Let an integral operator R^ admit the representation (3.3.11) 
where probability distribution fi satisfies the infinite-volume FK-DLR equa- 
tions (3.3.8). Then V g G G 

R [ ; ] = S^\g)R { ;\S^{g))- 1 . (3.3.13) 
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The proofs: the compactness and the tuned- 
action arguments 



The proof of Theorems 2.1 and 3.1 is based on a compactness argument (cf. 



[8], |9J). We want to note that this argument does not depend upon the 
dimensionality of the system. 



4.1. Proof of Theorems 2.1 and |3.1[ As in [8], [9], we first prove 
that, Vn>0, the sequences of RDMKs {Fjtff, N = n + 1, n + 2, . . .} and 

{-fftxw N = n+l,n+2, ...} are compact in the space C° (M xV " xM xV "). 
Applying Lemma 1.5 from [S] (this lemma goes back to [26J), we will obtain 
that the sequences of RDMs {R^ N } and {Rplf\ x cf N \} are compact in the 



trace- norm topology in T-L(n). This yields the statement of Theorem 2.1 



A straightforward consequence of the convergence will be that any limiting 
RDMK F^ admits the representation (3.3.4) where /i satisfies the infinite- 



volume FK-DLR equation, i.e., the assertion of Theorem 3.1 



To verify compactness of the RDMKs {F^} and 



we follow 



/3,N\x c (N) 

the same line as in [8], [9], i.e., employ the Ascoli-Arzela theorem. To this 
end, we need to check the properties of uniform boundedness and equiconti- 
nuity. For definiteness, we focus on the (slightly more complex) case of the 
sequence {F^ C{N) }. 

More precisely, to show uniform boundedness, we first use an upper bound 
for the number of vertices fc, on Vj \ Vj_i under the measure P; cf. [10], Eqn 
(4.1). Namely, V e G (0, 1), for P-a.a. T G Too 3 a constant C = C(T) such 
that 

,1/2+e 



ki < Ci(\n i\ 
This yields that 



2,3, 



(see HO] 

f; hj{i) < d(T) + c(t) ^(^) 1/2+e Q 



\ni 



(4.1.1) 



C{T)J\ (4.1.2) 



We use (4.1.1) and (4.1.2) to bound the quantity 
q(u>(n)\uj(N \ n)V * C (N)) := 



cxp 



- h(ut(n)\u{N\n) Vx c (iV)) 



l fitn (u(N\n) Vx c (iV)) 



(4.1.3) 
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cf. (3.3.7) for n' = n. Namely, (4.1.2) implies that, for P-a.a. T G Too, V 
n > and N > n, 

exp [-/3(U + C(T)J*F)#V n ] 



< exp 



h(u(n)\u(N\n) V x c (iV)) (4.1.4) 
<exp [/3(Z7 + C(T)J*F)ttV n ] 



for all path configurations a; G Wx(n),y(n) and loop configurations u)(N\n) G 
^Vjv\v„- Here (J V n — Y^=i ^ stands for the number of vertices in the set V n , 
cf. (4.1.1). 

The lower bound in (4.1.4) yields that 

E p , n (u(N\n) Vx c (iV)) > exp [-0(U + C(T) J*F)ft V n ] x (4.1.5) 
where 

PM-jAyT, E exp(-|n| 2 /2/3) (4.1.6) 

n=(n 1 ,...,n d )eZ d 

is the probability density of transition from x G M to x in time (3 in the 
Brownian motion on M. Next, (4.1.5) and the upper bound in (4.1.4) imply 
that 

q(u(n) \u(N \ n) V x c (iV)) < — ^ exp [2/3(Z7 + C(T) J*F)tf V B ] . (4.1.7) 

Substituting (4.1.7) in (3.3.7), we obtain that 

^|xc W (x(n),y(n)) < exp [2/3(17 + C(T) J*V)% V„] (4.1.8) 

which gives the desired uniform upper bound. 

To check equicontinuity, we analyze the derivatives 

v z« F /Sixc(7v)( x ( n )>y( n )) and v ?/(i) F ftjv|x-(jv)( x ( n )'y( n ))' * G v «- A s ain 

we use the representation (3.3.7) with n = n'. We need to differentiate the 
integral 

/ F f(n), y( n)( d ^)) ex P [-h(u(n)\u(N \ n) V x c (iV))]. (4.1.9) 
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For definiteness, consider one of the gradients V^j). It is convenient to 
represent the integral (4.1.8) in the form 



II tol-^T-k / <„),x(«)(<M«)) 



:(n),x(n) 



x exp | - /i[(w(n) +rj(n))\w(N\n) V x c (iV)] | . 
Here p] l/[ (x,y) denotes the transition probability density 



n={m,...,n d 



and p M has been determined in (4.1.6). 

Next, fj(n) = {r)(j),j G V n } is a collection of linear paths 

T 

vU, t) = -p{y{j) - X{j)), j G v n , 

and the component-wise addition in u>(n) + 77(71): 

w(n) + 77(71) = {w(j) + r/(j), j G V„} 

where w(j) + rj(j) : r G [0, j3] 1— >■ (w(j, r) + r/(j, r)) mod 1. 

It is now clear that there will be two contributions into 
V j(') f S|xc()V)( x W'y( n )) : one coming from 

the other from 

V„(i) exp { - /i (w(n) + 77(71)) |u> (TV \n) Vx c (iV) } 
= -V y(i) /i [(w(n) + 77(71)) I w(7V \ n) V x c (iV) 

x exp j - h [(w(ra) + 77(71)) I w(JV \ n) V x c (iV)] | . 

The uniform bound 

VvwAWUU)) <c(/3)g(o,+oo) 



(4.1.10) 



exp (-\ x -y + n\ 2 /2P) , (4.1.11) 



(4.1.12) 



(4.1.13) 



(4.1.14) 



20 



is straightforward. Next, we have the estimate 



(w(ra) + fj(n)) | w(iV \ n) V x c (iV) 

</3(#V n )[Z7 + C(T)J*F]. 



Together with (4.1.4) it implies that 



Vy(i) exp 



(w(n) + 77(71)) | w(iV \ n) V x c (iV) 



< {#V n )[U + C(T)J*V}exp [f3(U + C(T)J*V)$V n ] . 
The bounds (4.1.14) and (4.1.15) lead to a uniform bound upon 

^y(i) F p n) N \ x c ( N ) ( x M > y W ) 
EH □ 



(4.1.15) 



This completes the proof of Theorems |2.lj and 

4.2. Proof of Theorems 12.21 and 13.21 Theorem 12.21 follows from 
Theorem 3.2 therefore, we focus on the proof of Theorem 3.2 Eqn (3.3.13) 
follows from the property 



lim 



1, ge G, 



(4.2.1) 



uniformly in the path configurations u){n) = {co(j),j E V n } E W^f(n)y(n) an< ^ 

the loop configurations f2°(n') over V\V n '. Here, the functional q^\uJ(n) \ f2 c (n')) 
emerges from representation (3.3.11): 

q%\uJ(n)\n c (n'))(x(n)Mn)) = ^ 



&uj{n' \ n) 



x 



x(n),y(n) 



Pf (n)iy(n) (duf(n)) exp [-/i(5J(n) V «(n' \ n)| n c (n'))l , 

(4.2.2) 



and 



S{g)u{n) = {S{g)u(i),i EVn} 

where 5(g)w(z) : r e [0, /3] H- S(g)u(i,r). 

To check (4.2.1), we again follow the argument used in [H], [9] (which 
goes back to [16] and [1]; cf. also [5]). The backbone of the argument is 
the following inequality: V given a > 1, g E G and positive integer n, if n' 
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is large enough then, V u>(n) G y („) , x ( n )>y( n ) e M xVn and the loop 

configurations f2 c (n') over V \ V n ', 



aq%\ g ZJ(n)\n c (n>)) + a^VM")!*^"')) 



>2^(^(n)|fi c (n')). 



(4.2.3) 



The verification of Eqn (4.2.3) is based on a special construction related 
to a family of "tuned" actions g n n n u>(n' \n) on loop configurations u(n'\n); 
see Eqns (4.2.5), (4.2.6) below. A particular feature of the tuned action g n /\ n 
is that it "decays" to e, the unit element of G (which generates a "trivial" 
identity action), when we move the vertex of the tree T from V n towards 
V \ V n '- Formally, (4.2.3) is implied by the following estimate: V given n, 
uj(n) G W x ( n j y(n ), g G G and a G (1, oo), for any n! large enough, uj(n' \ n) 
and fT(n'), 



exp 



a 

+ -exp 
> exp 



h{gu(n)Vg nl> u(n'\n) \ n c (n')) 

-^g-^jVg^f^n) | n c (n')) 
-h(u(n) Vw(ra'\ra) | fl c (n')) 



(4.2.4) 



Indeed, (4.2.3) follows from (4.2.4) by integrating in du>(n' \n) and nor- 
malizing by S /3jn /(r2 c (n / )); cf. (4.2.2). Here it is important that the Jacobian 
of the map u>(n' \ n) i— >■ g n ,. \ n) is equal to 1. 

The rest of the argument concentrates on verifying (4.2.4). The tuned 
family g n /\ n consists of individual actions g^ n ' G G at vertices j G V n > \ V n : 



Sn'\n = {gj" ' 3 £ Vn'\V n }. 



(4.2.5) 



We use the representation (2.3.3) and identify the element g G G with a 
vector 6 = 9A<E IR '. Then the actions g^ n ^ G G correspond with multiples of 
the vector 9; cf. Eqn (4.2.6) below. It is convenient to fix a positive integer 
f > n and identify 

gf' ] with 9rf{n',k) (4.2.6) 
where = d(s , j) (recall, s is the root of T) and 



7(n', fc) 




k < r, 

f,n' — f), k > r. 



(4.2.7) 
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In turn, the function $(a, b) is determined by 

0(a, 6) = l(o < 0) + 1( ^ Q ~ / z(u)du, a,beR, (4.2.8) 



with the same functions Q(b) and z(u) as in [1] 



b 



(4.2.9) 



Q(6)= / *(u)du, 

Jo 1 

where 2; (it) = l(u < 2) + l(w > 2)— , 6 > 0. 

um u 

Next, g~/L is the collection of the inverse elements: 

g;V={gf r \.?eV n ,\Vn}. 

It will be convenient to use formulas (2.4.6)— (2.4.7) for g*- for j G V„, or 

even for j G V, as these formulas agree with the requirement that gj n ^ = g 

when j G V n and g^ n ^ = e for j G V \ V n ' . Accordingly, we will employ the 

notation g n , = {g\ j G V n 'j. 

Next, we use the invariance property (2.3.4). The Taylor formula for the 
function V G C 2 (IR 2 ) yields for j,f G V n : 



(4.2.10) 



v-(gfV^g?'W 

+^ (gf r ^0'),g?' r ^0")) - 2V(u;(j)Mf)) 

<C\9\ 2 \ 7 (n',j)->y(n',f)f i V. 

Here C G (0, 00) is a constant, the upper bound V is taken from (2.3.2), and 
we use the notation from (4.2.7). 

The bound (4.2.10) is crucial: this where the structure of the group action 
is exploited. It is based on the fact that the first-order terms in the expansion 
in the LHS of (4.2.10) cancel each other, due to the presence of elements g^ n ■* 

and ■* and their inverses, gj n ■* and gp ■* . This idea can be traced back 
to [IS] and [I]. 
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Further, the term \^{n' , j) — ^(n',j')\ 2 can be specified as 



| 7 (n',A;)- 7 (n',A;')| 2 =^ 



(4.2.11) 



0, if k, k' < f, or k, k! > n', 
[$(k — f,n' — r) 
— $(k' — f,n' — r)] 2 , 

if f < k,k' < n', 
•&(k — r,n' — r) 2 , 

if r < k < n', k' $]f, n'[, 
<d{k' — f,n' — f) 2 , 

if r < k' < n', k $]f, n'[ 



with notations k = d(j, s ), k' = d(j', s ). 

The convexity property of the function exp, together with Eqn (4.2.10), 
yield that, V a > 1, 



a 

+ -exp 



a - exp [ - h(g n , (uJ(n) V w(n' \ n)) |ft c (n')) 

" - h(g-}{uj(n) V w(n' \ n)) \fl c (n' 

-^( gri ,Hn)V^(n'\n))|fi c (n')) 



> a exp 



1 

2 

> a exp 



(4.2.12) 



-4(g n /(4)V W (n'\n))|n c (n'; 



4(«(n)Vw(n»|fiV)) e- c */ 2 



Here 



$ = $(n,g) = \9_\ 2 J WJ,J')) h(n,k) - 7 (n', A;')| 2 . (4.2.13) 

(i,i')ev„xv 

The series in (4.2.13) converges for P a.a. T, owing to condition (2.3.5) and 
estimate (4.1.16) below. 

The next observation is that 



$<3|£| 2 E l(k<k')j(d(j,j')) 

C?J')ev„,xV 



(4.2.14) 



x 



${k — r,n' — r) — <d{k' — r,n' — r) 
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where, by virtue of the triangle inequality, for all j,f : k < k! 



< d(k — r,n' — r) — $(k' — r,n' — r) 

,, z(k — f) 
< d(j,f)- 



(4.2.15) 



Q{n'-f)' 



Thus, 



$ < 



Q{n '- f)2 («0ev„,xv 

m 



J(d(jj'))d(j,f) 2 z(k-ry 



Owing to (2.3.5), it remains to bound the sum ^2 jeVn+ _ z(k — r) 2 . Note that 
u(\nu) 1 l 2+e z(u) < 1 when u G (wo(£),oo). Next, we use the bound (4.1.1) 
on the number of vertices in V n \ V n -i- Therefore, 



E z(k-r) 



E z(k — r) E z (k ~ r ) 



and 



j£V n '+r l<k<n'+r jeV k \V k ^ 

<C E #-r)<dQ(n'-f) 

l<k<n'+r 



C(T) 

$ < — „ ->■ oo, as n' ->■ oo. 



_ (4.2.16) 

Hence, given a > 1 for n' large enough, the term ae~ c */ 2 in the RHS of 
(4.2.12) becomes > 1. Consequently, 



exp 
a 

+ -exp 



- h (g n , (w(ra) V \ n)) | $~2 C (ra') 
> exp 



/i(a;(n) Vw(n'\n)|J]' 



V)) 



(4.2.17) 



Eqn (4.2.17) implies that 

gW(5i(n)|n c (n')) 



du?(n / \ n) 



(n'\n) 



X 



exp [ - h(p(n) V w(ra' \ n)|fi c (n))] 
S A n'(n c (n')) 
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(4.2.18) 



for any n and n 1 large enough obeys 



ri?>(g^(n)|n c (n')) + C{g-^{n)\n\n>)) 



(n), 



>2q%\uj(n)\n c (n>)) 



(4.2.19) 



uniformly in the boundary condition Q c (n'). Thus, (4.2.3) is established, 
which completes the proof of Theorem |3.2| □ 
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